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ABSTRACT 

We characterize the incidence of active galactic nuclei (AGNs) in 0.3 < z < 1 star- forming galaxies 
by applying multi-wavelength AGN diagnostics (X-ray, optical, mid-infrared, radio) to a sample of 
galaxies selected at 70- /im from the Far-Infrared Deep Extragalactic Legacy survey (FIDEL). Given 
the depth of FIDEL, we detect "normal" galaxies on the specific star formation rate (sSFR) sequence as 
well as starbursting systems with elevated sSFR. We find an overall high occurrence of AGN of 37±3%, 
more than twice as high as in previous studies of galaxies with comparable infrared luminosities and 
rcdshifts but in good agreement with the AGN fraction of nearby (0.05 < z < 0.1) galaxies of similar 
infrared luminosities. The more complete census of AGNs comes from using the recently developed 
Mass-Excitation (MEx) diagnostic diagram. This optical diagnostic is also sensitive to X-ray weak 
AGNs and X-ray absorbed AGNs, and reveals that absorbed active nuclei reside almost exclusively 
in infrared-luminous hosts. The fraction of galaxies hosting an AGN appears to be independent of 
sSFR and remains elevated both on the sSFR sequence and above. In contrast, the fraction of AGNs 
that are X-ray absorbed increases substantially with increasing sSFR, possibly due to an increased 
gas fraction and/or gas density in the host galaxies. 
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1. INTRODUCTION 

It has become increasingly clear that we need to 
reconcile the formation of stars in galaxies with the 
growth of the supermassive black holes at their cen- 
ters for a complete picture of galaxy evolution. The- 
oretically, active galactic nuclei (AGNs) are invoked as 
a means to quench star formation in galaxies thereby 
explaining the presence of massive red galaxies at re- 
cent epoch. Without this so-called AGN feedback, both 
semi-analytic models and cosmological si mulations tend 
to overproduce massive blue galaxies (e.g.. ICroton et all 
I2006t iGabor et all 1201 lh . Together with the tight rela- 
tion between supermassive black ho l e (SMBH) mass and 
bulge mas s (Magorrian et al.lll998t iFerrarese fc Merrittl 
l2000t iTremaine et al.l [2002?) , this leads to a picture of 
co-evolution bet ween galaxy stellar cont ent and central 
SMBH (but see lJahnke k, Macci6ll2011l for an alterna- 
tive interpretation). However, observational evidence of 
a physical connection between star formation and AGN 
remains sparse and mostly indirect, especially during the 
growth phase (pre-quenching). 

Different triggers of SMBH gro wth have been pro- 
posed: major galaxy mergers (ISanders et al.1 119881 : 
iDi Matteo et al.l 120051: iHopkins et all 120061) . large-scale 
disk instabilities (|Bournaud et al.l 1201 If ), and a secu- 
lar growth where SMBH growth is unrelated to the 
star formation r ate (SFR) of galaxies (jShao et alj|2010t 
ICisternas etafl 120111 : iMullanev et all l2012h . Some 
of these studies favor hybrid models where the 
high-luminosity end (quasar regime, with Lagn > 
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10 44 erg s 1 ) is dictated by major galaxy mergers while 
the l ower-luminosity AGNs follow a secular evolution 
(e.g.. IShao et all 120101: iLutz et all l2010t iRosario et all 
|2012|) . However, the role of AGNs and their underly- 
ing physical connection with host galaxies remain un- 
certain. Most studies of the host galaxy morpholo- 
gies of high-redshift (> 0.5) AGNs, which have tried 
to distinguish the relative importance of major merger 
and secular gr owth, have relied on X-ray selected sam- 
ples of AGN (iGabor et al.l [200l ICisternas et ail I201H: 
ISchawinski et al.ll201lUKocevski et al.ll2011[ ). The main 
drawback of this approach as a test of the major merger 
hypothesis is the potential insensitivity to a key phase 
of AGN growth - when the SMBH is expected to be 
deeply buried in the g as-rich center of the merging system 
(| Sanders et al.l [19881) . During that buried phase, X-ray 
emission from the AGN may be mostly absorbed by in- 
tervening material with high column densities and thus 
may be undetected in even the deepest X-ray surveys 
currently available. Thus, a real test of this connection 
and of the AGN content of actively star-forming galaxies 
requires the identification of both absorbed and unab- 
sorbed AGNs. 

Nebular emission lines tracing the narrow line re- 
gions (NLRs) surrounding AGNs are not subject to the 
same small-scale obscuration as X-rays owing to their 
much larger physical extent (100's of pc to a few kpc 
scales). Indeed, the [O in] A5007 luminosity combined 
with hard (2—10 keV) X-ray luminosity has been used 
as a Compton-thic kness parameter in order to infer the 
X-ray absorption ( Mulcha ev et all 11994 iBassani et al.l 
Il999t iHeckman et alj|2005t Uuneau et alj|2011i hereafter 
Jll; also see Appendix [SJ. AGNs identified at visible 
wavelengths tend to either have a direct view of the nu- 
clear region with low-level obscuration by dust in the 
host galaxy (type 1 AGNs, recognized by broad emis- 
sion lines) or a large obscuration of the nuclear region 
along the line-of-sight combined with a lesser obscura- 
tion of the narrow-line regions (type 2 AGNs, lacking 
broad emission lines 18 ). 

In extreme cases where both X-ray emission and opti- 
cal line emission are completely obscured, infrared light 
may reveal an AGN by showing reprocessed thermal 
emission from AGN-heated dust grains. Compared to 
the bulk of stellar light, the more energetic radiation 
field from an AGN heats the dust to higher temperatures, 
adding a hot dust compon ent with thermal emission at 
shorter MIR wavelengths (lElvis et all 11994 lLacv et al.l 
12001 iStern et all 120051: iMullanev et all 120 111) . This ex- 
tra hot dust alters the mid-IR spectral energy distribu- 
tion (SED) by producing an excess between the usual 
stellar bump at ~ 1.6 /im and the stellar- heated dust 
emission at > 10 ^m. Broad-band MIR observations 
such as those from Spitzer/IRAC can probe that feature 
(lLacy et all 12001 IStern~etaLl 120011: iDonlev et all 120071 
12012ft . Not all AGNs will exhibit MIR signatures be- 
cause of the required geometry and dust content. Also, 
intrinsically weaker AGN occurring alongside star for- 
mation tend not to dominate the IR SED of their host 

18 Additionally, some type 2 (narrow lines only) AGNs may 
be radiatively inef ficient and simply lack broad line regions (e.g. 
I Trump ct al. 2011) but these systems are intrinsically weaker and 
presumably have less effect on their host galaxies. 



galaxies ()Barmbv et all 120061 ) . 

Overall, it is clear that no single diagnostic can achieve 
a complete census of AGNs for all galaxies. As such, 
combining selections at multiple wavelengths may be the 
key to achieving both a better sampling and a better 
understanding of AGNs and their connection with star 
formation activity. 

In this paper, we characterize the incidence of AGN 
in star-forming galaxies at intermediate redshift (0.3 < 
z < 1). The multi- wavelength dataset and galaxy sam- 
ples are described in Sections [2] and [3l respectively. The 
AGN diagnostics are introduced in Sectional Section [5] 
contains the results regarding the AGN fraction among 
star- forming FIR-selected galaxies ( ij5.ip : the occurrence 
of active nuclei as a function of infrared luminosity ( £!5.3p 
and as a function of specific star formation rate (sSFRs; 
tj5.4[) : the absorption of AGN in relation to host galaxies 
( fl5.5P : and the influence of AGN emission on mid-to-far 
infrared color, probing the dust temperature on the warm 
side of the infrared SED ( tj5.6[) . The caveats are discussed 
in Section 16.11 followed by the possible connections be- 
tween AGN obscuration and host galaxies (Section 16. 2p . 
the properties of the AGNs selected via different wave- 
length regimes (Section l6.3p . and possible physical inter- 
pretations of the triggering mechanisms of AGN in Sec- 
tion [nUl Lastly, we summarize our findings in Section [7J 
We assume a flat cosmology wit h £L m = . 3, 17a = 0.7, 
and h = 0.7 throughout, and a IChabrien (|2003l) initial 
mass function (IMF) when deriving stellar masses and 
SFRs. 

2. MULTIWAVELENGTH DATASET 

The primary intermediate-redshift galaxy sample is 
based on observations from the Great Observatories Ori- 
gins Deep Survey 19 (GOODS) and the All-wavelength 
Extended Groth strip International Survey 20 (AEGIS), 
specifically in the GOODS-North and the Extended 
Groth Strip (EGS) fields. 

2.1. Infrared Photometry 

Both fields were observed during the Far- 
Infrared Deep Extragala ct ic Le gacy survey (FIDEL; 
Dicki nson fc FIDEL team! 120071) . yielding sensitive 
Spitzer observations at 24^m and 70/xm (respective 3a 
limiting fluxes of ^20/iJy and 2.5 mJy). MIPS photom- 
etry was obtained using a guided extraction method, 
where images at shorter wavelength, less subject to 
confusion, are used to build a set of prior positions to fit 
for the same sources at longer wavelengths. In this case, 
IRAC 3.6 /xm images were used to select priors which 
were then fit simultaneously at 24//m. Next, sources 
with a 24/zm detection were fit at 70jum. The data and 
method are described in more detail by Magnclli et all 
(120091 120111) . Using Monte Carlo simulations, these 
authors found that this method can deblend sources 
that are at least 0.5 x FWHM apart (~ 9" for MIPS 
70/im), and they quantified the systematic uncertainties 
using the difference between the ext racted flux of 
simul ated sources to the real input flux (Magnclli ct al. 
2009). The latest catalogs in GOODS (jMagnelli et~ 
120111 ) include two estimates of the uncertainties: the 

19 http://www.stsci.edu/science/goods/ 

20 http://aegis.ucolick.org/ 
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systematic uncertainties quantified from the Monte 
Carlo simulations, as well as local background noise 
uncertainties (from the residual maps). We use the 
maximum of the two values for each galaxy. 

Spitzer/IRAC photometry is available in all four 
channels (3.6, 4.5, 5.8 and 8.0/im; available through 
the Spitzer Science Center for GOODS-N and from 
iBarmbv et all [20081 for ECS). In what follows, IRAC 
photometry is used to calculate rest-frame K-band mag- 
nitudes, which are in turn used to estimate stellar masses 
when the latter are not available through SED fitting. 
In addition, IRAC photometry is used for mid-infrared 
color-color AGN diagnostics (£ I4.3|) . 



2.2. Optical Spectra 

As in the work by lJuneau et al.l (|2011l hereafter 
Jll), the GOODS-N optical spectra are drawn from 
the T eam Keck Redshift Survey 21 (TKRS IWirth et all 
120041 ). whereas the EGS spectra come from the DEEP2 
Galaxy Redshift Survey (hereaft er DEEP2; iDavis et"al"1 
120031 12007t iNewman etalll2012l) . However, we augment 
the sample with spectra obtained durin g the DEEP3 
camp aign includin g galaxies in GOOD S-N (jCooper et al.1 
l201lh and EGS ([Cooper et al.1 l20ll Cooper et al, in 
prep.). 

All sets of observations were obtained with 
Keck/DEIMO S and reduced with the sam e pipeline 22 
(|Cooper et al.l l2012t INewman etldl 12012ft . However, 
their spectral resolution and spectral range differ due to 
the use of different gratings (600 line mm -1 for TKRS 
and DEEP3, and 1200 line mm" 1 for DEEP2). The 
TKRS and DEEP3 resolution is 4A FWHM over the 
wavelength range 5500-9800A, wher eas DEEP2 spectra 
have a resolution of 2A FWHM with a wavelength 
coverage of 6500-9100A. 

Emission line fluxes are measured as described in Jll, 
by fitting a Gaussian curve to each line individually and 
integrating the flux over ±2.5<r, where a is the Gaussian 
width (=FWHM/2.355). Through this paper, we are 
only concerned with the H/3/[0 in] A5007 flux ratio and 
with the [O in] A5007 luminosity. The H/3 fluxes are cor- 
rected for underlying Balmcr absorption cither by using 
a population synthesis mo del fit to the continuum (with 
iBruzual fc Charlotl (|2003h models when the median S/N 
per pixel is greater than 3), or by applying the median 
equivalent width of 2.8 ±0.9 A found in Jll to correct for 
Balmcr absorption for spectra with low S/N. 

2.3. X-ray Data 

We take advantage of the fact that the Chan- 
dra X^^_coverage is ve ry deep: 2 Ms in GOODS- 
N (lAlexander eta D 12003ft . 200 ks (|Laird et al.l 2009; 
iNandra et al.ll2005l) . and 800 ks in the EG S (Laird et 
al., in prep.; reduced in a similar fashion as lLaird et al.1 
l2009ft . In the EGS, we favor the 800 ks X-ray fluxes 
when available, and use the wider area 200 ks observa- 
tions otherwise. The shallower EGS X-ray data (200 ks) 
are sufficient to detect moderate luminosity X-ray AGNs 
{L2-wkeV > 10 42 erg s _1 ) out to z ~ 1. Furthermore, 

21 http: // tkserver . keck, hawaii .edu / tksurvey / 

22 http://deep.berkeley.edu/spec2d/, developed by the DEEP2 
team at the University of California-Berkeley 



we can detect fainter X-ray galaxies (mostly starbursts) 
in GOODS-N and in the 800 ks portion of the EGS. 

X-ray observations are primarily used for AGN identifi- 
cation (Section |4Tj. Furthermore, they are used to infer 
X-ray absorption by comparing the hard X-ray luminos- 
ity (£2-iofceV, calculated as in Jll) to the [O in] A5007 
narrow- line region luminosity (L^q m i), a more isotropic 
tracer of the intrinsic AGN luminosity than X-rays. For 
that purpose, we derive X-ray upper limits for [O ni]- 
sel ected AGNs. Up p er lim its were derived as described 
by lAlexander et~aLl (|2003h in the GOODS-N field, and 
estimated from the sensitivity maps in the EGS as de- 
scribed by Jll. 

2.4. AGN Bolometric Luminosity 

Bolometric AGN luminosities are estimated using 
Z/boi(AGN)= 10 3 xL[ O in i. The bolometric conver- 
sion factor to L2--10 koV has been calculated to range 
from 10—70 depe nding on Eddington ratio and /o r 
X-ray luminosity ([Vasudevan &: Fabianl l2007l l2009ft . 
Here we assume Lboi(AGN) = 25 x L 2 -io keV and an- 
other factor of 40 to convert from [O in] to unab- 
sorbed X-ray (2— lOke V) luminosity (the value found by 
iHeckman et al.l ([2005) for unobscured AGNs). 

The AGN bolometric luminosities can be converted to 
a BH mass accre tion rate with the followin g equation 
(Equation 1 from lAlexander fc Hickoxl (|2012ft ): 



M BH [M Q yr- i ]=0.15( — 



^bol 



10 45 ergs" 1 



(1) 



where e is the efficiency of conversion of mass into en- 
ergy. We assume a value e = 0. 1 in our calculations 
([Marconi et al.ll20p: lMerlonill200l . 

2.5. Infrared Luminosity and Star Formation Rates 

Total infrared luminosities (8 < A < 1000/im) were ob- 
tained by fitting the 70 /im flux d ensity and redshift with 
templates from ICharv fc Elbad ([200ll hereafter CE01). 
Even when used with only data at much shorter wave- 
lengths (e.g., MIPS 24/im) than the peak of typical in- 
frared SEDs (^60-100 /im), these templates have been 
shown to give consistent results with direct integration 
over the longer wavelength Herschel bands that probe 
the peak of the F IR SED in galaxies (out to z ~ 1.5; 
lElbaz et al.ll20ToT) . 

Because we are interested in the IR luminosity origi- 
nating from star formation only, observations at 24 /im 
or shorter wavelengths were not included in the fit. Rel- 
ative to 70 /im, they risk to have a larger contamination 
from AGN-heated dust when galaxy nuclei are active (see 
Section [6. ip . Using a contaminated 24 /im measurement 
would lead to an overestimate of the IR luminosity due 
to stellar heated dust and of the SFR. It is possible that 
even observations at 70 /im have an AGN contribution 
but the effect will be less important than at shorter wave- 
lengths. 

Total infrared lumi nosities we r e conv erted to SFRs us- 
ing the relation from iKennicuttl (11998ft with a addi tional 
conversion fact or of 0.66 from a ISalpeterl fl955) to a 
iChabrierl ((2003ft IMF. 

3. GALAXY SAMPLES 
3.1. Intermediate- Redshift Samples 
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We start from a parent sample of spectroscopically 
selected galaxies at z > 0.3 with a spectral coverage 
that includes the expected wavelengths of the H/3 and 
[O III] A5007 emission lines. The resulting redshift range 
is approximately 23 0.3 < z < 1. We will refer to this 
parent sample of 9435 galaxies as the Inter/ All sample 
(where Inter stands for intermediate redshift). 

The Inter/FIR subsample is further defined to have a 
MIPS 70/im detection with S/N> 3, yielding a sample 
of 270 galaxies at 0.3 < z < 1 (58 in GOODS-N; 212 
in EGS). Galaxies in the Inter/FIR sample all have a 
robust 24 /im detection (S/N> 8) given the much greater 
sensitivity of the 24 /im observations relative to the 70 /im 
ones. The Inter/FIR sample includes 94% of MIPS 70 
galaxies with a spectroscopic redshift 0.3 < z < 0.8. 
This percentage decreases to 82% over the redshift range 
0.3 < z < 1 because the limited spectral range of the 
DEEP2 setup prevents the coverage of the H/3 and [O in] 
lines beyond z > 0.8. 

Among the Inter/FIR sample, 145 (54%) galaxies have 
valid emission line measurements with a 3tr detection for 
both the H/3 and [O III] A5007 lines, while 62 (23%) 
galaxies have one detection and one upper limit. This 
results in 207 (77%) galaxies with a constraint on their 
[O in] A5007/H/3 flux ratios. The remaining galaxies have 
S/N< 3 for both emission lines. 

Given the sensitivity of the FIDEL survey at 70/im, 
the majority (71%) of the 70/Ltm-selected galaxies at 
0.3 < z < 1 are luminous IR galaxies (LIRGs; with 
Lir(8 - 1000/xm) > 10 11 L Q ), as shown in Figure [Q 
These objects are responsible for a large fraction of the 
global SFR density at z ~ 1 (|Le Floc'h et al.l 120051: 
iMagnelli et al.l I2009D and extend down to the typical 
star- forming galaxies on the main-sequence of sta r forma- 
tion (jNoeske et al.ll2007t lElbaz et al.ll2007tl2011l) . Thus, 
they comprise not only highly starbursting systems, 
which are thought to only contribute ~ 10% of the global 
star formation rate dens ity (SFRD) (|Rodighiero et al.l 
1201 It ISargent et aT]|2012D . but also more representative 
galaxies that contribute to the remaining 90% of the 
SFRD. 

3.2. Low-Redshift Comparison Sample 

We build a z ~ comparison FIR-selected sam- 
ple from the Sloan D igital Sky Survey (SDSS) DR7 
(|Abazaiian et al.l 120091) based on the detection in the 
IRAS 60^011 or AK ARI 90/im bands as described by 
I Hwang et alj (|2010al fbl). Following the procedure in Jll, 
we further select galaxies to be at redshift 0.05 < z < 0.1 
and to have emission line detections with S /N> 3 for all 
features used in the spectral classification from Jll (H/3, 
[O in] A5007, Ha, [N n] A6584, and [S n] AA6717, 6731) 
yielding a sample of 4034 galaxies, which we call the 
Low/FIR sample. 24 

The emission line fluxes were obtained from the Value 

23 In detail, there is some variation in the wavelength range of 
individ ual s pectra due to the instrumental set-up (as described in 
Scction l2.2jl and due to the exact position of the slits on the mask. 

24 The S/N requirement was relaxed for galaxies with SDSS 
flag spclass = 3 in order to retrieve broad-line AGNs. Among 64 
spclass = 3 objects (1.5% of Low/FIR sample). 40 have S/N> 3 in 
the catalog and the remaining 24 either had a lower S/N or were 
not measured in the MPA/JHU catalog due to the breadth of the 
lines (> 700 km s _1 ). 
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Fig. 1. — Total infrared luminosity as a function of redshift for 
70/im detected galaxies in GOODS-N and the EGS (all symbols). 
Galaxies with spectral coverage of both the H/3 and [O III] A5007 
emission lines are shown with filled squares (Inter/FIR sample), 
and the galaxies lacking coverage are shown with open squares. 
These features shift beyond the wavelength range of DEEP2 spec- 
tra by z ~ 0.8, explaining the smaller fraction of galaxies with 
available emission lines at 0.8 < z < 1. The X-ray classification 
indicates robust X-ray AGNs (open circles) and X-ray ambiguous 
cases (open star symbols) with faint and soft X-ray emission like 
X-ray starbursts. We show the 3<r sensitivity limit of the FIDEL 
survey at 70/im with the shaded area and the 3<r limit and 80% 
completeness limit of the 24/im data (dashed and dot-dashed lines, 
respectively). All the limits were derived with the CE01 templates. 

Added Catalogs developed by the Max-Planck Institute 
for Astronomy (Garching) and John Hopkins Univer- 
sity (MPA/JHU) 25 . The methodology fo r measurements 
of em ission line fluxes is described by iTremonti et al.l 
( 2004). We adopt the AGN classification based on the 
BPT [N il]- and [S li]-diagrams as described by Jll, but 
using the classification scheme from lKewlev et all (j2006[ ) 
does not alter our results. 

4. MULTIWAVELENGTH AGN DIAGNOSTICS 

4.1. X-ray Diagnostic 

The X-ray identification of AGNs is based on two crite- 
ria: (1) hard X-ray luminosity L 2 -iokeV > 10 42 erg s _1 , 
or (2) hardness ratio 26 HR > —0.1, corresponding to an 
effective photon index T < 1. These conditions are simi - 
lar to the X-ray criteria employed by I Bauer et al.l (|2004[ ) 
based on observed trends that nearby galaxies without 
AGN do not exceed a lumin osity of 3 x 10 42 erg s _1 in 
the 0.5 - 8 keV band (e.g., iFabbianol l"l989D . and that 
X-ray sources of stellar origin are characterized by com- 
paratively soft X-ray emission with pho ton indices r > 1 
(IColbert et al.ll2004T also see Figure 2 of lAlexander et all 
2005). 

The 2— lOkeV rest-frame luminosities were computed 
using ^-correction as described by Jll but not cor- 
rected for intrinsic absorption within the galaxies. Sev- 
eral other studies also adopt such a luminosity threshold 

25 http:/ /www. mpa-garching.mpg.de/SDSS/DR7/ 

26 Hardness Ratio = (H - S)/(H + S), where H and S are the 
X-ray count rates in the hard (2 — 8 keV) and soft (0.5 — 2 keV) 
bands, respectively. 
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TABLE 1 

Definition of Galaxy and AGN Samples 



Sample Number Redshift Criteria 

Low-Redshift Comparison Galaxy Sample 

Low/FIR 4,034 0.05 < 2 < 0.1 S/N>3 for H/3, [O in] A5007, Ho, [N n] A6584 and [S n] AA6717, 6731 emission lines 

and detection in IRAS 60 /im or AKARl 90 fim 

Intermediate-Redshift Galaxy Samples 

Inter/All 9,435 0.3 < z < 1 TKRS and DEEP2/3 spectra that cover the wavelengths of H/3 and [O III] A5007 

(from parent samples selected with iJ[AB]< 24.3 and 24.1, respectively) 
Inter/FIR 270 0.3 < z < 1 Inter/All galaxies with a valid detection with MIPS 70 fim (S/N>3) 

Intermediate-Redshift AGN Selections 

X-AGN 19 0.3 < z < 1 L 2 _i fc e v > 10 42 ergs" 1 and/or HR>-0.1 

IR-AGN 15 0.3 < z < 1 AGN region of IRAC two-color or single-color diagram 

MEx-AGN 92 0.3 < z < 1 MEx-AGN probabilities 

Radio-AGN 5 0.3 < z < 1 Radio-excess (q 70 <1.6) 

Intermediate-Redshift AGN Samples 

X-Unabsorbed 17 0.3 < z < 1 X-ray identified AGNs a 

X-Absorbed 64 0.3 < z < 1 X-ray unidentified AGNs b but intrinsically luminous 

[log(L [0 „ l] [ergs- 1 ])>40.4 or IR-AGN or radio-AGN] c 

Weak 18 0.3 < z < 1 Fail X-AGN, IR-AGN, radio-AGN selections and with weak [O in] A5007 

[log(L[o in] [erg s- 1 ]) < 40.4] c 



"Excluding the two most extremely absorbed X-AGNs with log(Lx/£[o rn])< 0-3. 
'Including the two systems with X-AGN identification but with \og(Lx / Liq In ])< 0.3. 

c Using the same conversion factors as in i|5.2l log(L[o ml) = 40.4 corresponds to log(Z/2_io kev) = 42 and to log(Li, ;) = 43.4 
with luminosities in units of ergs -1 . 



when selec ting samples of AGNs from deep X-ray sur- 
veys (e.g., ICisternas et alj|2011t ISchawinski et aJJl2011t 
iMullanev et al.l r2012f ) . The addition of a hardness cri- 
terion in the current study allows us to include weaker 
and/or more absorbed AGN, although only one among 
19 X-ray AGN has L 2 - Wk ev < 10 42 egs s _1 . 

X-ray detected sources that do not fulfill the AGN cri- 
teria are considered "X-ray ambiguous", as the X-ray 
emission could have a stellar origin rather than being due 
to AGN activity. However, a faint and soft X-ray emis- 
sion can also be present in extremely absorbed AGNs if, 
for instance, there is a concurrent X-ray starburst or if 
the X-ray spectrum is reflection-dominated rather than 
transmission-dominated. Nearby Compton-thick AGNs 
are observed to span a broad range of spectral indices 
(see Figure 11 of I Juneau et alJl2011h . 

4.2. MEx Diagnostic Diagram 

The second AGN identification method is the MEx di- 
agnostic diagram introduced by Jll. The MEx diagram 
plots stellar mass versus the [O in] A5007/H/3 flux ratio 
(Figure Ek). Empirical dividing lines indicate spectrally 
distinct regions of the MEx diagram: optically-identified 
AGN tend to lie above the lines, purely star-forming 
galaxies below the lines and the MEx-intermediate re- 
gion (between the lines) contains a sharp transition from 
star-forming to AGN galaxies as indicated by the strong 
gradient from P(AGN)=0% to P(AGN)=100%. While 
we use the dividing lines for discussion and visualization, 
the analysis relies more importantly on the AGN prob- 
ability values 27 . Those probabilities are determined by 
using a set of priors from the SDSS as described by Jll 
and correspond, for a given observed galaxy, to the frac- 
tion of AGNs of any type (Seyfert 2, composite, LINER) 
relative to the total number of galaxies located within 
the measured uncertainties on the [O III] A5007/H/3— M* 

27 See https://sites.google.com/site/agndiagnostics/ for soft- 
ware to compute the MEx AGN classification. 



plane. 

The main strengths of the MEx diagram are its ap- 
plicability to higher redshift relative to traditional BPT 
diagrams 28 and the use of a probabilistic approach to 
assess the reliability of a given spectral class. In the 
remainder of this work, we use the AGN probability 
P(AGN) to count the number of AGN (e.g., a galaxy 
with P(AGN) = 30% is counted as 0.3 AGN instead of 
1 AGN) in cases where an AGN is identified only with the 
MEx method. The AGN fraction is defined as follows: 

1 N 

AGN fraction = — ^ P(AGN). t (2) 

i=l 

where the AGN probability P(AGN) varies from to 
1 and is summed over the number of galaxies N that 
belong to the subsample of interest. This approach re- 
duces the risk of contamination by star-forming galax- 
ies into the AGN samples but may lead to undestimat- 
ing the total AGN fraction as composite galaxi es, which 
mostly harbor AGNs (e.g. iTrouille et al.ll201ll) . tend to 
have P(AGN)<100%. The initial values of P(AGN) are 
obtained from the MEx diagnostic (Pmex(A-GN)). How- 
ever, AGNs identified from an alternative method (X- 
ray or IR colors) are set to have P(AGN) = 1. The 
68.3% confidence intervals on fractions are determined 
with Bayesian binomial statistics following the fo rmal- 
ism and IDL 29 implementation of ICameronl (|2011h . 

The MEx diagnostic diagram successfully detects 
AGNs that are missed in even the deepest X-ray sur- 
veys. Jll showed that at least some of the X-ray un- 
detected galaxies with a probability > 30% of hosting 
an AGN contain X-ray absorbed AGN activity. Indeed, 
X-ray stacks revealed detections in both the soft and 
hard bands of the Chandra X-ray observations with a flat 

28 BPT diagrams require Ha and other fairly red emission lines 
([N n] A6584, [S n] AA6717, 6731, or [O i] A6300), which redshift 
out of the visible window at z > 0.4. 

29 Interactive Data Language 
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spectral slope indicating X-ray absorption. Although the 
MEx AGN classification may not always hold on an in- 
dividual galaxy basis, it was v alidated overall to z = 1 
for average galaxy populations (puneau et al.ll201l"t ) , and 
preliminary studies at slightly higher redshifts suggest 
that the MEx method remains valid on a yerage (z ~ 1.5; 
Trump et al., submitted). Furthermore. iMullanev et all 
(2012) found that the mean stellar mass of X-ray se- 
lected AGN hosts did not evolve between 0.5 < z < 3. 
This trend suggests that the diagnostic power of the MEx 
diagram may hold at those higher redshifts for the mod- 
erate to luminous AGNs (L^oi > 10 43 4 ergs -1 ). Never- 
theless, we limit the analysis to z < 1 in this work, where 
the MEx diagnostic can be applied with a high level of 
confidence. 

4.3. IRAC Single- Color and Two-Color Diagnostics 

Infrared AGN identi fication method s rely on probing 
AGN- heated dust. lLacv et al.l f|2004h and IStern et al.l 
(|2005h reported that combinations of IRAC colors se- 
lect both obscured and unobscured luminous AGN. 
These studies were based on fairly shallow IRAC data. 
iBarmbv et "all (|2006l) found that the Lacy and Stern MIR 
diagnostics, when applied to more sensitive IRAC data, 
miss a significant fraction of X-ray selected AGNs and 
suffer from more contamination from non-AGN galaxies. 
The latter is caused by the deeper observations probing 
to higher redshift (z > 1) and making the IRAC colors of 
starburst and normal star-forming galaxies appear redder 
and moving into the AGN selection box ()Donlev et al.l 
l200l . 

Here, we adopt the IRA C color-color diagram devel- 
oped by IStern et al.l (|2005f ) because it was shown to suf- 
fer from less star-forming galaxy contamination than the 
Lacy diagram when applied to hig her-redshift and t o 
more sensitive IRAC observations ()Donlev et alj 12007). 
The diagram is shown in Figure [2f b) , where AGNs nomi- 
nally lie in the upper boxed region. We show some tracks 
of expected mid-infrared colors for various galaxy tem- 
plates. We do not expect strong contamination by star- 
forming galaxies at z < 1, but some well-known AGNs 
are located outside of the AGN region when their colors 
are redshifted. 

Because IRAC color-color diagrams do not include in- 
formation about the redshift of the galaxies, there are 
degeneracies in observed colors from overlapping red- 
shift tracks of diff erent galaxy populations (see, e.g., 
iDonlev et all 12007). To remove such degeneracies, we 
use the IRAC [3. 6] -[8.0] color as a function of redshift. 
As shown in Figure EIc), infrared-selected star- forming 
galaxies occupy a distinct region from nearby AGN 30 , es- 
pecially at z > 0.5. The outermost star- forming galaxy 
template is for IRAS 12112+0305, which we use to sep- 
arate the IRAC star-forming galaxies (below the cyan 
curve) from the IRAC AGN candidates (above the cyan 
curve). The change in IRAC [3.6]-[8.0] color with redshift 
of the infrared selected star-forming galaxies is due in 
part to the aromatic features at rest-frame wavelengths 
6.2 — 8.6 /im gradually shifting outside of the [8.0] band 
(between < z < 0.4) and the stellar bump at rest-frame 

30 The AGN templates were constructed by gathering photome- 
try from the NASA Extragalactic Database and interpolating with 
templates (R. Chary, private communication, 2008). 



~ 1.6 /xm entering the [3.6] band at 0.7 < z < 1.2. When 
AGN-heated dust dominates the mid-IR emission, the 
near power-law shape of the continuum causes a flatter 
[3.6]- [8.0] change of color with redshift, as exemplified by 
Mrk 231 and UGC 5101. Lastly, dust- free s tellar popu- 
lations templates (Br uzual k Ch ariot 2003) set a lower 
envelope to the observations. When galaxies are nearly 
or completely dust-free, the IRAC color corresponds to 
that of stellar photospheres. 

4.4. Radio AGN Diagnostic 

Radio-excess AGNs can be identified by an excess of 
radio emission relat i ve to the FIR-radio correlation (e.g. 
Rov k Norriail997t IDonlev et ah! 12001 iDel Moro etaLl 
2012D . which otherwise is very tight for star-forming 
galaxies and many radio-quiet AGNs (see review by 
ICondonl 11992). Radio AGN identification was found to 
reveal a different AGN population with limited over- 
lap with the X-ray and IRAC color AGN diagnostics 
(|Hickox et al.l f2009) . However, radio AGNs were also 
found to reside in different host galaxies, characterized 
with, on average, higher stellar masses, stronger cluster- 
ing, and passive stell ar activity compar ed to the hosts 
of X-ray or IR AGNs (|Hickox et alj|2009D . The typically 
more passive galaxies expected to host radio AGNs are 
unlikely to be detected in the MIPS 70 /im observations 
used to draw our sample. Therefore, we expect that this 
category of AGNs will not be important for the sample 
selected here. 

In the GOODS-N field, there are deep Very Large 
Array (VLA) 1.4 GHz observations re aching 3.9 iiJy 
beam rms with a beam size of ~1."7 JM 

orrison et al] 

2010). T hese data were u s ed to identify radio-excess 
AGNs by IDel Moro et al] (|2012l ) for counterparts of 
24 /im sources which have S/N>3 in the 1.4 GHz obser- 
vation (5(1.4 GHz) > 12/jJy). Infrared SED fitting was 
used to compute FIR emission integrated over rest-frame 
~ 40—120 /im in order to apply a standard FIR-radio cor- 
relation 31 . Only one galaxy in our GO ODS-N subsam- 
ple is selected as a radio-excess source bv lDel Moro et all 
(|2012t ) . This galaxy is also identified as AGN from all the 
other methods (X-ray AGN, IR AGN, and MEx-AGN) 
and therefore does not add a new AGN population. 

In the EGS, VL A observations were presented by 
(jlvison et alj [2007) along with a catalog of sources 
with 5(1.4 GHz) > 50 (i3y (~5cr) and a ~3."8 FWHM 
synthesized beam. For the 41 radio sources with a 
match in the Inter/FIR sample, we computed q^o = 
log(57o fim/SiA GHz) as a proxy for the FIR-radio cor- 
relation. The median of (770=1. 9 is consistent to values 
previously report ed in the same redshift range (1.9±0.2 
and 1.9±0.1 from lFraver et alj|2006t iBourne et al]|201lL 
respectively) and, as exp ected, slightly lowe r than k- 
corrected q*m values (see I Sargent et al] 120101 and the 
compilation in their Table 5). Adopting q-j^ < 1.6 to 
define radio-excess AGNs yields four AGN candidates 32 . 

31 This was shown by the authors to be more sensitive 
than using MIPS 24 fim as a m onochromatic p roxy of the FIR 
with the selection criterion of IDonlev et alj 1 20051) : 924 = 

log(S 2 4 fim/S 1A GHz) < 0. 

32 This cut corres ponds to a >3 a excess according to the work of 
IBourne et al.l 1 20111), and to a 1.5 a excess according to the work 
of IFraver et alj 120061 ). The significance of the excess may thus 
depend on the number of galaxies and the sample selection. 
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Fig. 2. — AGN diagnostic diagrams: (a) The MEx introduced in Jll. The empirical curves (solid lines) split the diagrams into galaxy 
spectral types as labeled, with a MEx-intermediate region between the MEx-AGN and MEx-SF classes. Logarithmic contou rs show the 
numb er density of the SDSS galaxy sample used to define AGN probabilities (color bar), (b) IRAC color-color diagram from Ste rn et al.l 
(2005). The AGN selection box (thick solid line) and the power-law locus (thin solid line) are marked in gray. IRAC colors are computed 
for various templates as a function of redshift from z = 0.3 (open diamo n ds) to 2 = 1. Namely, we display the expected colors of nearby 
star-forming LIRGs and ULIRGs from the templates from Rickc ct al. (2009, in yellow), nearby AGNs in light blue, dust-free stellar 
population models using BC03 templates in purple, and we highlight a nearby Compton-thick double AGN, NGC 6240, with a thick black 
line, (c) IRAC color between observed channels 1 and 4 as a function of redshift. The cyan curve, which corresponds to the redshift track of 
IRAS 12112+0 305, is the star-form ing galaxy template showing the most extreme IRAC colors. The remainder of the star-forming galaxy 
templates from [Rickc ct al. (2009) (yellow tracks) span a fairly restricted range in color. Galaxies above that range are suspected to host 
AGN activity. This picture is supported by the projected colors of a few nearby AGNs [black lines: Mrk 231 (solid); UGC 5101 (dotted 
line); and Mrk 273 (dashed line); but also see Compton-thick AGN NGC 6240 (thick solid)]. Dust-free stellar population models from 
BC03 define the lower envelope (purple tracks). 



Two of them are identified as MEx-AGNs (but unde- 
tected in X-rays and classified as star-forming on the 
IR diagnostics), while the other two are not identified 
as AGN by any method. One of the unide ntified ra- 
dio AG Ns has optical signatures according to lYan et aT] 
(|2011h . who assumed a more relaxed contraint on emis- 
sion line detection (S/N>2 instead of S/N>3 as in the 
current work). 

Combining both GOODS-N and EGS fields gives a to- 
tal of five radio-excess AGN candidates, three of which 
are also identified with alternative methods. The remain- 
ing two radio AGNs make only 0.7% of the Inter/FIR 
sample, and f.9% of the AGN subsample. Therefore, 
their inclusion or exclusion does not alter the results pre- 
sented in this Paper, but we include them in the analysis 
whenever the plotted quantities are available. Namely, 
the two radio- AGNs that are otherwise unidentified will 
be put in the X-ray absorbed category (Section 15. 2[) . 

5. RESULTS 
5.1. AGN Identification 

The Inter/FIR sample is shown on the MEx diagram 
in Figure Eta). Most (16/19=84%) of the X-ray AGNs 
lie in the MEx-AGN and MEx-intermediate regions, as 
expected. In addition, there are 78 AGNs identified from 
the MEx diagnostic that were missed in the X-ray clas- 
sification, including 14 X-ray ambiguous sources. 

In Figure [3] (b) , we compare IRAC color-color iden- 
tified AGN with the X-ray and optical MEx classifica- 
tions. The X-ray AGNs display a similar distribution 
as the optically-selected AGN on the IRAC color-color 
plane. A few of the IRAC-AGN candidates which are 
undetected in X-rays are classified as AGNs with the op- 
tical diagnostic. Combining both the X-ray and MEx 
classifications selects 90% (9/10) of the galaxies within 
the IRAC selection box. 



The X-ray, MEx and IRAC color-color AGN selec- 
tion methods are next shown on the IRAC single-color 
diagnostic described in Section 14.31 and displayed in 
Figure EJc). Unsurprisingly, the FIR-selected galaxies 
mostly follow the tracks of the nearby infrared galaxies 
from Figure [21 Galaxies hosting X-ray AGNs (black cir- 
cles) span the whole range of IRAC colors covered by the 
Inter/FIR sample (from non-AGN and up to the most 
extreme AGN-like color). The X-ray AGNs that are in- 
distinguishable from star-forming galaxies are likely sys- 
tems where star formation dominates at infrared wave- 
lengths even in the pr esence of nuclear activity (e.g., 
iRujopakarn et al.ll2011l ). We note additional AGN candi- 
dates that were not selected by the X-ray method (open 
black circles) nor by the IRAC color-color diagram (Panel 
b). In the remainder of this paper, IRAC-selected AGN 
(or simply IR-AGN) refer to AGN identified from either 
the color-color or the single-color versus redshift method. 

5.2. Summary of AGN Selection and AGN Properties 

We illustrate the overlap and complementarity be- 
tween the AGN classification methods in Figure |4j The 
overlap between the MEx method and the IR and X- 
ray methods is greater for the more luminous AGNs, 
as expected because they are easier to detect with any 
given tracer explored here. X-ray selection suffers from 
ambiguity between faint AGNs and starburst galaxies 
at low X-ray luminosities (£2-10 keV < 10 42 ergs _1 , 
or L b oi(AGN)< 10 43 4 ergs^ 1 ) and the IRAC selec- 
tion is mostly sensitive to intrinsi cally luminous and 
dusty AGNs (Lw > 10 43 erg s" 1 ; iHickox et all 120091 : 
iDonlev etT al. 2012). Although optical emission lines can 
be subject to attenuation by dust in the host galaxy, 
they remain the most sensitive AGN tracer and probe 
to much lower accretion rates due to the detection of 
relatively faint emission lines in optical spectra (in this 
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Fig. 3. — AGN diagnostic diagrams introduced in Figure [2] (a) The MEx AGN diagnostic diagram. The empirical curves (solid lines) 
split the diagrams into galaxy spectral types as labeled, with a MEx-intermediate region between the MEx-AGN and MEx-SF classes. 
Contours show the SDSS low- 2 sample (evenly spaced logarithmically), and the 70/^m sample is superimposed with larger symbols keyed to 
the P(AGN) (color bar) and to the presence of X-ray AGN (circles) or ambiguous X-ray source (possibly X-ray starbursts, star symbols), 
(b) IRAC color-color diagram from [Stern et af] f 20051 ) . The colored symbols are keyed to the MEx-AGN probabilities (color bar) whereas 
the gray plus symbols mark galaxies lacking a valid MEx classification due to low S/N [O III] A5007 and H/3 lines. X-ray identified AGNs 
(X-ray ambiguous sources) are further marked with a circle (star symbol). Most IRAC AGN candidates in our sample are also recognized 
as such from the X-ray and/or MEx selection methods, (c) IRAC color as a function of redshift. Same plotting symbols as the previous 
panel. The lower lines are the BC03 tracks for dust-free stellar populations whereas the top gray lines is the template of IRAS 12112, used 
as an AGN-star formation dividing line on this diagram. 



work, down to fluxes of ~ 7 x 10 -18 ergs -1 cm -2 , cor- 
responding to Lboi(AGN)^ 10 42 5 ergs _1 and accretion 
rates of 5 x 10~ 4 M© yr _1 ). The greater sensitivity likely 
explains the additional AGN populations that are not 
detected with alternative methods. However, there may 
be cases where the narrow line region emission lines are 
extinguished by surrounding dust. It is important to 
combine different AGN methods to take advantage of 
their complementarity and achieve a more complete cen- 
sus. Furthermore, the combination of multiwavelength 
information can distinguish between different AGN sub- 
populations (or AGN accretion phases). 

The different AGN tracers are combined in order to 
divide the AGN population into three categories that are 
physically meaningful: 

1. X-unabsorbed AGN: X-ray identified AGNs 
except for two systems that have the high- 
est absorption (Njj > 10 24 cm -2 , inferred from 
log(L x /L [0 m] ) < 0.3; Appendix© 33 . 

2. X-absorbed AGN: AGNs that are unidentified 
in X-rays - plus the two most absorbed systems 
excluded from the X-unabsorbed category - but 
which are inferred to be intrinsically luminous due 
to their identification with an IR or radio method 
or due to having Lr [n l A5007 > 10 40 4 ergs -1 , cor- 
responding to Lx > 10 42 ergs -1 . These AGNs 
should thus be detectable with the X-ray method 
in absence of absorption. 

3. Weak AGN: AGNs which are not identified with 
X-rays nor IRAC colors and which have a faint 
[O ill] A5007 luminosity (L [Q m] < 10 40 - 4 ergs -1 ), 
implying that they would not be recognized as 



MEx 




Radio 



X-ray 



33 The remainder may contain moderately absorbed AGN as we 
do not measure the absorbing column density from X-rays. The 
latter would be highly uncertain due to low number of counts for 
most objects. 



Fig. 4. — Venn diagram of the AGN subsamples: MEx-AGN (in 
green), radio-excess AGN (in purple), IR-AGN (in orange) and X- 
AGN (in black). The shaded areas are proportional to the number 
of galaxies in each category and overlap region. The MEx-AGN 
sample is further split into bins of AGN bolometric luminosity 
(Lf, z) as labeled. The luminosity cut at Lj, D ; > 10 43 4 ergs - iwould 
only remove one AGN from the IRAC AGN subsample, correspond- 
ing to 7% (1/15) of the area shown in the figure. 

AGNs with the X-ray definition even if they were 
X-ray detected. 

These AGN subsamples are mutually exclusive, and the 
subsample sizes and definitions are summarized in Ta- 
ble [TJ They share interesting similarities and differences 
in their host galaxy and accretion rate properties, and 
may correspond to different phases in a typical SMBH 
growth cycle (Section 153)) . 

In what follows, we compute the distributions of AGN 
and star-forming (SF) galaxy subsamples by adding the 
probabilities of each galaxy falling in the category of in- 
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terest. The number of galaxies per bin is the sum of (1- 
P(AGN)) for the SF subsample and the sum of P(AGN) 
for the AGN subsamples. When calculating numbers, 
P(AGN) is converted as a number from to 1 rather 
than a percentage. Similarly, the probability-weighted 
means are obtained as follows: 

N 

^I,x P{AGN), 
(X) = ^ (3) 

where X is the quantity of interest (A/*, £rom]j or 
L[o m]/Lm), and the sums are over all the N galaxies in 
a given subsample. P(AGN) is replaced with l-P(AGN) 
for the SF sample. This methodology takes full advan- 
tage of the probabilistic approach of the MEx diagram. 

As can be seen in Figure [SJ X-unabsorbed AGNs are 
only found in the most massive hosts (~ 10 10 ' 6 ~ 12 M©) 
whereas X-absorbed and weak AGNs span a broader 
range of stellar masses (~ 10 10 ~ 12 M©) and have a simi- 
lar average stellar mass as the non-AGN galaxies in the 
Inter /FIR sample (~ 10 10 - 6 M Q ), which is ~ 0.5 dex 
lower than the average stellar mass of the X-unabsorbed 
sample. The apparent lack of X-unabsorbed AGNs in 
lower mass hosts may be due to a combination of AGN 
evolutionary phase and selection bias (discussed in Sec- 
tion MM- 

Absorbed and X-unabsorbed AGNs have similar AGN 
luminosities on average (Lboi(AGN)^ 10 44 ergs _1 ; in- 
ferred from [O in]; Figure [5j and are systematically 
higher (by 1 dex on average) than the weak AGNs and 
the SF galaxies. The physical meaning of Lboi(AGN) 
does not apply to the SF subsample (unless the galax- 
ies contain unidentified AGNs). However, the similarity 
of the [O ill] luminosities between weak AGNs and SF 
galaxies (similar means, and relatively similar distribu- 
tions) implies that we are able to distinguish between 
these two populations based on an additional factor and 
that the [O in] luminosity alone is not sufficient to dis- 
criminate between AGNs and star-forming galaxies. This 
may be because [O in] may contain a contribution from 
low-metallicity gas ionized by star formation. This con- 
trib ution is expected to be negligible in metal-rich galax- 
ies (jKauffmann et al1l2003fl , but we test this possible con- 
tamination to [O in] by examining the [O in]-to-IR lu- 
minosity ratios. If the [O III] emission traced the SFR 
instead of the AGN luminosity, and assuming that Ljr 
traces the SFR, we would expect similar ratios between 
star forming galaxies and the false AGN candidates. In 
the case of genuine BH growth with small or no stellar 
contribution to [O III], the [O ni]-to-IR ratios can be 
converted as black hole accretion rate (BHAR)-to-SFR 
ratios. Absorbed and X-unabsorbed AGNs have similar 
BHAR-to-SFR ratios, and are systematically higher (by 
0.5-0.6 dex) than the weak and SF subsamples indicat- 
ing an [O in] excess emission (Figure [5|). The later two 
subsamples exhibit similar [O ili]-to-IR distributions to 
one another, suggesting that the weak AGNs are not only 
weak relative to the absorbed and X-unabsorbed subsam- 
ples, but also relative to the SFR of their hosts. Some 
of the weak AGN may represent a distinct AGN regime 
with possibly radiatively inefficient accretion (|Ptak et al.l 



119981: lHoll2008t iTrump et aTll201lD rather than a simple 
extension of the same accretion pr ocess down to lower 
luminosity (but also see lMaozll2007f ). 

5.3. AGN Fraction in Intermediate- Redshift Infrared 
Galaxies 

In order to obtain the global AGN fraction, we com- 
bine all of the AGN identification methods, namely: the 
MEx diagnostic diagram, the X-ray classification, and 
the IRAC single-color and two-color diagrams. Adopting 
the procedure from Section 14.21 and Eq. 14.21 the global 
AGN fraction in 0.3 < z < 1 FIR-selected galaxies is 
37%±3% (99/270). 

The AGN fraction increases steadily with Lm and 
both the increase and the magnitude of the AGN fraction 
are consistent with the values obtained for the low/FIR 
comparison sample (Figure [6j . The similarity between 
the AGN fractions of the Low/FIR and Inter/FIR sam- 
ples is consistent with no or very mild evolution in the 
triggering of AGN as a function of IR luminosity between 
z ~ 1 and z ~ 0. 

In contrast, the AGN fraction found for the Inter/FIR 
sample is significantly higher than th at of similarly se- 
lected galaxies in the COSMOS field (jKartaltepe et al.1 
l2010al) . Those authors combined X-ray, IRAC power- 
law, and radio AGN identification methods (also see 
ISvmeonidis et ai1l20 10) but lacked a robust optical emis- 
sion line diagnostic 34 . Our higher AGN fractions are 
mostly a result of a more complete AGN selection due 
to the inclusion of X-ray faint and absorbed AGNs from 
the MEx diagram. The latter account for most of the 
difference between the samples (right panel of Figure [5]) . 

There are more X-ray absorbed AGN than X-ray un- 
absorbed AGN at all IR luminosities but the distribu- 
tions of the infrared luminosities of these two AGN cat- 
egories are similar. The intrinsically weaker AGN tend 
to reside in less infrared-luminous hosts. The dearth of 
intrinsically weak AGN in the most luminous galaxies 
may be due to increased dilution by star formation and 
not necessarily caused by a difference in the triggering 
mechanisms. If we missed AGN due to star formation 
dilution, our AGN census would be more incomplete at 
the IR-luminous end. However, we already observe a 
larger AGN fraction in more IR-luminous systems, which 
is opposite to the effects of dilution or dust obscuration. 
Thus, the observed increase in AGN fraction with Lir is 
robust despite this potential bias. 

5.4. AGN on the SFR-M* Sequence 

The high incidence of AGN in infrared-luminous hosts 
may hint at a link between black hole growth and SFR. 
The SFR is known to evolve with redshift and to corre- 
late strongly with stellar mass at a given redshift. The 
correlation between the SFR a nd stellar mass has been 
reported in several studies (e.g.. lBrinchmann et alj|2004l : 
iNoeske et aLll2007t lElbaz let al.H2007| ). Assuming a lin- 
ear relation between SFR and stellar mass at a given 
redshift 35 , the evolution of the normalization has been 

34 ISvmeonidis et al.l <201Ch used a simple [O III] A5007/H/3 > 3 
threshold, which only selects AGNs with the most extreme line 
ratios. This selection criterion would yield only 10 out of 92 (11%) 
MEx-AGNs. 

35 The index a of the relation SFR oc M" is consistent with unity 
within the uncertainties and within the spread of values published 



10 



Juneau et al. 



log L bol (AGN) [erg s" 1 ] 
42 43 44 45 



log M BH /SFR 
-4 -3 



30 



& 

w 



20 



10 



SF 

X-unabsorbed 

weak 



■ ■■■■■■■■ L 



~Vr .f77~i . . i 




10 11 12 

log stellar mass [M s ] 



39 40 41 42 

log 4om] [ er g s "'] 




-5 -4 -3 

lo g L [011I]/ L 1R 



Fig. 5. — Distribution of host galaxy stellar masses (left), [O in] luminosities (center; corresponding AGN bolometric luminosities on 
top axis), and [O lll]-to-IR luminosity ratios (right; corresponding BHAR to SFR ratios on top axis). Those properties are shown for 
star-forming galaxies (blue dotted line) and the three AGN subsamples : X-unaborbed (black), X-absorbed (orange) and intrinsically weak 
(green). The physical meaning of Lf J0 i(AGN) and BHAR does not apply to the SF subsample (unless some of those galaxies contain 
unidentified AGNs). The number of galaxies is computed with 5~l (1 — P(AGN)) for the SF subsample and with J])P(AGN) for the AGN 
subsamples (as a reminder, P(AGN)= Pm Bx (AGN) for MEx-only AGNs and P(AGN)= 1 for X-AGNs and IR-AGNs). The probability 
weighted means (Eg. I5.2t are indicated with vertical lines. 
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Fig. 6. — Left. AGN fraction as a function of Lm in intcrmediate-redshift FIR-selected galaxies (filled histogram). The AGN 
fraction rises with in creasing Lju in agreement with the trend observed for low-redshift (0.05 < z < 0.1) FIR-selected SDSS galaxies from 
IHwan g et alJ (|2010al . shown in blue). For comparison, we also show values obtained for the 70/^m selected sample from Kartaltcpc ct al. 
J2010al) , spanning similar luminosity and redshift ranges (filled black symbols). Right. Similar to the left-hand panel except that the 
AGN content is broken down into X-ray unabsorbed, X-ray absorbed, and intrinsically weak systems, as labeled. Interestingly, the overall 
trend of the AGN fraction against Ljr resem bles that of the ga l axy m erger fraction (dashed line) an d th e combined irregular and merger 
morphology fraction (dotted line) observed by Kartaltcpc ct al. (2012) at z ~ 1 (discussed in Section 16.411 . Horizontal error bars show the 
bin width and v ertical error bars show the 68.3% confidence intervals for fractions, determined with Bayesian binomial statistics following 
the formalism of Cameron (2011). 
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parameterized bv lElbaz et al.l f|2011h as: 

sSFR sequence (z) [Gyr- 1 ] = 26 x t^ lc (4) 

where t cosm i C is the cosmic time in Gyr at the redshift z of 
interest. Here, we study the specific star formation rate, 
sSFR, of the AGN hosts relative to the remainder of the 
inter/FIR sample. The X-ray unabsorbed AGNs tend to 
reside in galaxies on or slightly below the sSFR s e quenc e, 
in agreement with results from Mullanc y et al.l ((2012)), 
while the X-ray absorbed and/or X-ray weak AGNs are 
hosted by galaxies with higher sSFRs (Figure [7]) . 

The overall AGN fraction appears to be independent of 
both the sSFR and Alog(sSFR), the relative offset from 
the redshift-dependent sSFR sequence: 

A \og(sSFR) = \og(sSFR) - \og(sSFR sequence (z)) (5) 

(Figure [5]) . Contrary to the near constancy of the AGN 
fraction with respect to the sSFR, we find a strongly 
increasing fraction of X-ray absorbed AGNs among 
the moderately-luminous AGN population as sSFR (or 
Alog(sSFR)) increases. Combined, these results suggest 
that the AGN obscuration depends strongly on the sS- 
FRs of the host galaxies, but the total incidence of AGNs 
(probed by the AGN fraction) does not. 

5.5. X-ray Absorbed AGN 

As shown in Section 15.31 and summarized in Table [TJ 
X-ray absorbed AGNs are more common in actively star- 
forming (FIR-selected) galaxies (24%) than are X-ray un- 
absorbed AGNs (7%). In what follows, we define the 
absorbed AGN fraction as the number of X-absorbed 
AGN over the combined number of X-absorbed and X- 
unabsorbed AGN, i.e., without taking into account the 
intrinsically weak systems. The behavior of black hole 
growth and AGN absorption is investigated as a func- 
tion as host galaxies properties. Figure |9] shows that the 
occurrence of AGN increases with infrared luminosity, 
tracing SFR, and is fairly flat with stellar mass down 
to 1O 1O M0. The uncertainty in the lowest mass bin is 
too large to constrain the behavior of the AGN fraction 
at masses < 1O 1O M0. The AGN fraction against the 
sSFR is consistent with a constant value around ^35- 
40% on and outside the sSFR sequence. Lastly, the AGN 
fraction increases with warmer 24-to-70/Ltm colors as ex- 
pected from an AGN contribution to dust heating (Sec- 
tion IQj). 

The behaviour of the absorbed AGN fraction differs 
from that of the occurrence of AGN in galaxies. The 
fraction of AGNs that are absorbed is fairly constant 
with IR luminosity with perhaps a mild increase at the 
highest L/fl, while it decreases with increasing stellar 
mass, and increases sharply with increasing sSFR. These 
trends suggest that AGN X-ray absorption may not be 
as sensitive to the total gas mass (somewhat traced by 
Liu) as to the gas fraction and/or gas density, traced 
by sSFR. The absorbed fraction decreases with warmer 
24- to- 70 yum colors, suggesting a smaller relative contri- 
bution of star formation in less absorbed AGN. Lastly, 
the absorbed fraction decreases with increasing bolomet- 
ric luminosity of the AGN, in agreement with results re- 
ported previously (e.g.lUeda et a l. 2003; La F ranca et al.l 
120051 iTrump et all I2009T) . iGilli et al.l (|200l proposed 
an analytical description of the ab sorbed fract i on ad - 
justed to represent the data points of lUeda et al.l (|2003|) : 



lLa Franca et al.l (|2005l) with their model m2, shown in 
the bottom right panel of Figure [5] 

5.6. Effect of AGN on Mid-to-Far Infrared Color 

The observed 24-to-70/mi flux ratios of the Inter/FIR 
galaxy sample span a broad range of values (Figure ITTjf , 
exte nding beyond t he ra nge of the IR SED templates 
from lDale fc Heloul (|2002l . hereafter DH02), especially at 
lower redshift 36 . We use the DH02 template that splits 
the sample in comparable numbers on either side as a 
redshift-dependent dividing line between galaxies with 
high and low S24/S70 galaxies. The chosen template has 
parameter a = 1.85, where a is the index of the power- 
law relating the mass of dust heated to the intensity of 
the interstellar radiation field responsible for the heating, 
see Equati on 1 of DH02. Accord ing to the method de- 
scribed by Marcill ac et al.l (|2006h . this SED corresponds 
to a total infrared luminosity of 10 11 Lq. X-ray identi- 
fied AGNs tend to reside above the line and reach the 
largest values of S^/SVo- 

With Figure [TTJ we study the mid-to-far infrared color 
of the AGN hosts (including X-ray, IR, and MEx AGNs) 
with respect to the rest of the star-forming galaxies in 
our sample. Galaxies with a low S24/SV0 (red triangles) 
occupy predominantly the star-forming and intermedi- 
ate regions of the MEx diagram (Panel a), and the star- 
forming region of the IRAC two-color diagram (Panel 
b). Galaxies with a high S24/SV0 (blue squares) are dis- 
tributed across the full range of parameter space with 
the most extreme cases (filled black squares) almost ex- 
clusively classified as AGN on the MEx diagram, and 
in or near the AGN region on the IRAC two-color dia- 
gram. Galaxy mid-to-far-infrared colors are more evenly 
distributed on the IRAC single-color diagram (Panel c) 
but there is an excess of the most extreme infrared col- 
ors (filled black squares) in the AGN region. Overall, 
galaxies with the highest S24/ SVo values are much more 
likely to be identified as AGN with any (or all) of the 
AGN tracers, especially X-rays and optical lines. This 
result indicates that the active nuclei within these galax- 
ies significantly contribute to the heating the dust whose 
emission is measured by the 24 and 70 /im passbands at 
these redshifts. 

The AGN-heated contribution to dust emission is ex- 
pected to peak at shorter wavelengths (~ 10 — 20 /im) 
compared to dust heated by star formation processes, 
and hence produce an enhanced mi d- to far-IR ratio . 
The average AGN SED temp late of lElvis et all (Il994f) 
suggests this feature (see also iMullanev et al.ll2011l for 
the intrins ic AGN IR SEP of m ore moderate luminos- 
ity AGN). IVeilleux et~aTl (|2002D found that the spectra 
of nearby ULIRGs with warmer IRAS 25-to-60 /jm color 
(S25/S60 > 0.2) are quasar-like whereas the spectral fea- 
tures of cooler ULIRGs are similar to LINER or purely 
star-forming galaxies. Similarly, we observe that the 
galaxies with the highest S24/SV0 values were selected 
as X-ray AGN (Figure HO]) , MEx- AGN (Figure EJi) and 
mostly also as IRAC AGN (Figure ITTb. c). The expecta- 
tions appear to hold well for these few extreme systems, 

36 The DH02 models were constructed using IRAS sources, and 
so are deficient in galaxies whose FIR is dominated by cold dust 
such as those sele cted in longer wavebands with, e.g., Herschel 
llSmith et al.ll20ll) . 
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Fig. 7. — Left. Specific star formation rate as a function of redshift. The inter/FIR galaxy sample spans the sSFR sequence, shown 
with a gray shaded area, and beyond. X-ray unabsorbed AGN (black asterisks) mostly lie on or below the sSFR sequence whereas X-ray 
absorbed AGN (orange circles) and intrinsically weak AGN (green squares) tend to have higher sSFRs. The latter are shown only for objects 
with an AGN probability percentage > 30% but the right hand panel conveys the information on the individual P(AGN) percentages. 
Right. Same as left panel, except that the symbols are color-coded according to P(AGN) when available (otherwise shown with black plus 
symbols). The shape of symbols indicate the [O III] luminosity (filled circles if > 10 40 ' 4 ergs -1 ; filled squares otherwise). X-ray unabsorbed 
AGNs are furthermore marked with large open circles. Contours indicate the distribution of the Low/FIR sample on the sSFR-z plane. 
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Fig. 8. — Left. Fraction of galaxies hosting active nuclei as a function of their specific star formation rate (red diamonds). Right. 
Fraction of galaxies hosting active nuclei as a function of their distance from the specific star formation rate sequence (red diamonds). In 
both panels, the AGN categories are distinguished with the shades of gray as labeled. The X-ray absorbed fraction among the moderately- 
luminous AGN population (e.g., excluding the weak AGN category), is shown with the open circles and dashed purple line. Horizontal 
error bars show the bin width and v ertical err o r bars show the 68.3% confidence intervals for fractions, determined with Bayesian binomial 
statistics following the formalism of ICamero 3 goU). 



but whether there is also a trend for the more normal, 
less extreme AGNs is less clear. 

Figure rT2] presents another diagnostic by combining 
the [O in] A5007/IL3 and S 2 aI 'S 70 ratios. There are 
202 galaxies in our sample with all required measure- 
ments, including galaxies with an upper limit for cither 
[O in] A5007 or H/3 (but not both). The bulk of the 
sample tends to occupy the bottom left part of the plot 
with several outliers toward the upper right, which are 
the more extreme and easily identified AGN (according 
to at least one of the X-ray and the MEx classification 
schemes) . 

We define an empirical curve to divide the region dom- 



inated by X-ray AGNs from the rest of the sample 

0.6 



log([0 m]/H/3) 



log(S 24 /S 7 n) + 0.28 



1.2 (6) 



A dividing line at fixed [O m]/H/3 ratio (log([0 m]/H/3)= 
-0.15) separates the optically selected AGNs (green to red 
colored symbols) from the star-forming galaxies (purple 
to blue colored symbols) . Given the fact that the galax- 
ies between the lines are not detected in hard X-rays and 
given their range of [O iii]/H/3 values, they may be com- 
posite SF/AGN galaxies and possibly LINERs (though 
we e xpect very few LIN ERs given the FIR selection; see, 
e.g.. lYuan et al.ll2010t ). Interestingly, the X-ray unde- 
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Fig. 9. — AGN fraction (red diamonds) and fraction of the AGNs that are X-ray absorbed (filled black circles) as a function of the 
following parameters: infrared luminosity, specific SFR, distance from the sSFR sequen ce (AsSFR) , stella r mass, mid-to-far infrared color, 
and AGN bolometric luminosity. In the last panel, the X-ray obscured fraction from [Gilli et al.l {2003) is shown with the shaded area, 
where we applied the same bolometric correction factor to the X-ray luminosity as in the rest of this work (1.4 dex). Horizontal error bars 
show the bin width and v ertical error bars show the 68.3% confidence intervals for fractions, determined with Bayesian binomial statistics 
following the formalism of Cameron 1 20111 ). 

tected AGN candidates and star-forming galaxies appear 
to have a very similar distribution in their 24-to-70/im 
colors. These systems may have a weak AGN whose 
contribution to the dust emission is small, or a more lu- 
minous AGN concurrent with an elevated SFR and thus 
also resulting in a small relative contribution. Alterna- 
tively, the IR SED shape might have a more complicated 
dependence on AGN heating owing to the geometry of 
the system. More generally, this result implies that fac- 
tors other than the presence of AGN should be taken into 
account in order to explain the spread in S^/SVo values. 

This analy sis is analogous to that presented by 

iKewlev et al.l (|2Q0lD . where the authors derived a mix- 
ing line going from low values of IRAS 25/60 flux 
ratio and [O m]/H/3 ratio toward higher values as the 
AGN fraction increases. Furthermore, they found AGN 
fractions of 73% and 77% in galaxies with k^S^/S^o) 
> -0.6 and -0.5, respectively. Although probing some- 
what different rest-frame wavelengths, our results are 
consistent with the same general picture, as all galax- 
ies with log(524/S , 7o) > -0.6 satisfy our AGN selection 
criteria. The results of increasing incidence of AGN 
in systems wi t h war mer mid-to-far infrared co lors from 
IKewlev et all (|2001h and IVeilleux et ail (|2002t) were de- 
rived for nearby galaxies. Figure [12] is the higher-redshift 
analog. 

The overall trend between the presence of AGN and the 
mid-to-far-IR color S24/ '>SVo is shown in Figure [TBI where 
the AGN fractions are calculated using the same method- 
ology as in Section 15731 We find an increasing AGN frac- 
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Fig. 10. — Observed 24-to-70 fira flux density ratio as a func- 
tion of redshift for galaxies in GOODS-N and EGS. The observed 
values are compared directly to redshifted model templates from 
I Dale fc Heloul 1)200 2'). applicable to star-forming galaxies (gray 
lines). The black line shows the DH02 template chosen to split 
our sample between high S'24,/S?o (blue squares) and low S^/SVo 
ratios (red triangles), with the most extreme S24/S70 highlighted 
with filled black squares. X-ray AGNs and X-ray ambiguous 
sources arc respectively marked with open circles and star sym- 
bols. For visual comparison, we include redshifted tracks for three 
nearby AGNs: NGC 1068 (Seyfert 2, in green), Mrk 231 (Seyfert 
1, in orange), and UGC 5101 (obscured or buried AGN, in purple). 
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Fig. 11. — AGN diagnostic diagrams introduced in Figure[2] Symbols are keyed to the mid-to-far-infrared color S24/S70 (red triangles 
for low values; open blue squares for high values; filled black squares for very high values) and to the X-ray classification (open circles 
for X-AGN and open star symbols for X-ambiguous). (a) The MEx AGN diagnostic diagram. The empirical curves (solid lines) split the 
diagrams into galaxy spectral types as labeled, with a MEx-intermediate region between the MEx-AGN and MEx-SF classes. Contours 
show the SDSS low-z sample (evenly spaced logarithmically). Objects with low values of S24/S70 tend to be centrally distributed on the 
MEx diagram while the location of the high-S^/SVo galaxies extends to both extremes of the [O m]/H/3 range. Galaxies with very high 
S2A/S10 are exclusively on the AGN tail of the distributio n toward high values of [O m]/H/3. (b) IRAC color-color di agram with the 
AGN region (thick solid line) defined by (Stern ct al. 2005), and the power locus (thin solid line) by ([Donley ct al. 2007). Galaxies with 
the highest S24/S70 are located within or nearby the boundaries of the AGN region, as expected if there were a 24/im excess caused by 
AGN-heated dust emission, (c) IRAC color as a function of redshift. The AGN candidates lie above the dark gray dividing line and span 
a broad range of S24/S70 values although 72% (13/18) have high or very high ratios (open and filled squares). The light gray tracks at the 
bottom define the lower envelope expected from BC03 stellar population templates with no dust emission. 
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Fig. 12. — Combined optical/IR diagram showing the 
[O III] A5007/H/3 emission-line flux ratio as a function of a mid- 
to far-IR color (24-to-70/xm). The 70^tm galaxy sample was classi- 
fied based on the MEx diagram (color keyed to P(AGN) as shown 
on color bar) as well as from their X-ray properties (star symbols 
for X-ray starbursts/ambiguous systems, and black circles for X- 
ray AGNs). The lines represent an empirical division between the 
most extreme AGNs (top right), the more common AGNs (middle) 
and star- forming galaxies (bottom left). 

tion with warmer S24/ '£70 values. While the warm color 
side appears to be driven by X-ray detected AGNs, re- 
moving these objects from the sample reduces the linear 
Pearson coefficient only from 0.94 to 0.88. Overall, we 
thus recover the expected trend from dust heating by 
AGN in the sample of FIR-selected galaxies used in this 
work. 

6. DISCUSSION 



Fig. 13. — AGN fraction as a function of mid- to far-IR color 
(24-to-70/xm). The FIDEL 70fim galaxy sample was classified as 
AGN based on all of the AGN diagnostics used in this paper: X- 
rays, MEx diagram, and IRAC c olor s. The split into AGN sub- 
categories is described in Section 15.31 and is identical to the right 
panel of Figure [6] Horizontal error bars show the bin width and 
vertical error bars show the 68.3% confidence intervals for frac- 
tions, deter mined with Baye sian binomial statistics following the 
formalism of Cameron (2011). 

6.1. Caveats 

There are three main caveats in this analysis although 
none of them is expected to significantly alter our re- 
sults. The first is that total infrared luminosities, which 
we interpret as measures of SFR, contain mixed contri- 
butions from dust heated by both stars and AGN. While 
longer wavelengths (70/im) are less subject to AGN con- 
tamination than, say, 24/im, they can still have non- 
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zero contribution from AGN-heated dust emission. Re- 
cent Herschel studies have shown that the FIR colors 
of galaxies with an AGN do not differ from those of 
purely star-formin g galaxies (|Hatziminaoglou et al.l2010t 
lElbaz et al.ll2010T ). indicating that on average AGN con- 
tamination does not dominate at A (observed) > 70 /zm. 
We mitigate against this effect by not using the shorter 
wavelength data (A < 24 /im) when calculating Lm for 
the FIR-selected sample (Section [53]). 

The second caveat concerns the overall completeness of 
the AGN selection methods used in this work. As men- 
tioned in Section 15. 3[ some objects have an unknown 
classification because of the low S/N of their emission 
lines. As a consequence, the AGN fraction that we de- 
rive for the overall FIR-selected sample (i.e., with and 
without emission line detections) may be underestimated 
as we consider these unclassified galaxies to be purely 
star-forming unless they fulfill the X-ray or IR AGN 
selection criteria. Taking into account that some of 
the unknown class are AGN would strengthen our con- 
clusion that AGN are ubiquitous among FIR-selected 
galaxies given that we already find a large AGN frac- 
tion. We computed a maximum AGN fraction assuming 
that all of the unknown class galaxies host an AGN: it 
would rise from 37% to 62%, but this scenario is un- 
likely given that the unknown galaxies do not fulfill the 
emission line criteria and that emission-line galaxies were 
shown to have a higher inciden ce of AGN than their non 
emission-line counterparts (e.g. lYan et al.ll2006T ). We test 
this trend with our own sample by comparing the X- 
ray and IRAC AGN fractions for the global FIR sample 
and the emission-line FIR sample. These methods sug- 
gest a 2x higher AGN fraction in emission-line galaxies. 
There are U±%% (18/145) X-AGN and IRAC- AGN in 
the emission-line subsample, but only 6±2% (7/125) in 
the remaining subsample, resulting in a global fraction 
of 9% (25/270) when considering all FIR-selected galax- 
ies regardless of the detection of [O in] and H/3 emission 
lines. Among objects with AGN classifications, it is also 
possible that the MEx diagnostic diagram misses AGNs 
preferentially in low-mass galaxies, especially if there is 
also ongoing star formation diluting AGN emission sig- 
natures. Such AGNs are fairly rare in the nearby uni- 
verse but would not be accounted for if they were more 
common at higher redshift. 

The third caveat concerns the use of [O in] A5007 as 
an indicator of bolometric AGN luminosity. On the one 
hand, [O in] emission could include a contribution from 
star formation yielding an overestimate of the total AGN 
power. T his is a more seriou s conc ern for low-metallicity 
galaxies. iKauffmann et all ()2003l ) found that in high- 
metallicity galaxies, the [O in] A5007 line is the least 
contaminated of the strong emission lines measured in 
optical spectra, with a flux fraction of 7% from star 
formation and the remaining 93% from AGN induced 
emission. The stellar m ass-metallicity (M* — Z) relation 
(jTremonti et al.l l2004f ) furthermore implies that metal 
poor galaxies have small stellar masses. Given that all 
the FIR-selected galaxies are fairly massive, they likely 
have a fairly elevated metallicity and their [O III] lines 
are likely dominated by AGN emission. However, the 
mass-me tallicity relation has been shown to evolve with 
redshift (|Savaglio e t al. 2005]) and/or SFRs of the galax- 



ies (|Mannucci et al.ll2010f) in the sense that the galaxies 
in the current sample could have a slightly lower metal- 
licity than galaxies of like stellar masses in the nearby 
universe. We thus compared L[ Q ml to Lju, tracing the 
SFR, in order to check whether the MEx- AGNs present 
[O in] excess due to their elevated SFRs rather than true 
AGN excitation. Both X-ray absorbed and unabsorbed 
AGNs show the same [O ni]-to-IR excess with respect to 
the subsamples of Weak AGNs and star-forming galaxies 
(Figure 0) . 

On the other hand, [O in] A5007 emission may be af- 
fected by dust obscuration. Even though we use [O in] 
as a more isotropic tracer than 2—10 keV, it is not i n 
fact perfectly isotropic (e.g. lDiamond-Stanic et al.ll2009h . 
Dust obscuration would have the opposite effect, and 
would lead to an underestimate of the intrinsic AGN lu- 
minosity when not applying a correction for extinction. 
We conservatively choose not to correct for dust obscu- 
ration, which means that the X-ray absorption may be 
underestimated and similarly, the number of absorbed 
AGN may be a lower limit. Given that we already find 
a high incidence of absorbed AGN among FIR-selected 
galaxies, our result would only be strengthened if we in- 
deed underestimated the number of absorbed AGNs. 

6.2. Link between AGN Obscuration and Host Galaxies 

A high absorbed fraction among AGNs at high red- 
shift may be expected because in addition to small-scale 
(torus) absorption, th ere were important gas reservoirs 
in ga laxies at z > 1 (jDaddi et al.1 120101 : iTacconi et al.l 
2010) that can potentially contribute to absorbing X- 
rays. Models of gas-rich galaxies predict high column 
densities (reaching the Compton-thick regime with Nh > 
10 24 cm -2 ) along several lines-of-sight in both isolated 
gas-rich unstable di sks dBournaud et al.l 1201 lh and in 
major mergers (e.g., iHopkins et al.ll2006j) . Furthermore, 
analyses of the cosmic X-ray background infer an im- 
portant population o f X-ray absorb ed AGNs at z < 1 
i|Comastri et~allll995t iMerlonil [200l . Therefore, a high 
fraction of hidden AGNs should be expected at increas- 
ing redshift. 

Is the absorption occurring on small scales or galaxy 
scales or both? In the former case, one would expect that 
AGN absorption would not correlate with galaxy-scale 
properties 37 . Interestingly, we found that the fraction of 
the AGNs that are absorbed increases sharply with the 
sSFRs (Figure [9]). Before discussing the physical impli- 
cations, we further examine X-ray absorption in relation 
to the infrared luminosities of the host galaxies with an 
extension to fainter host galaxy samples. 

We start from AGNs that have both X-ray and [O ill] 
detections in the Inter/ All sample (dubbed the X+[0 in] 
sample), and divide this population according to their 
detections in the MIR only (FIDEL 24 /im) or in the 
FIR (FIDEL 70 fim). The 24/zm comparison sample 
has a much fainter (~ 10 x) Ljr sensitivity limit than 
the 70 /xm sample (Figure [1} . The remainder of the X- 
ray+[0 III] AGN sample lacks FIDEL 24 /im detections 
and thus probes yet lower IR luminosities. These three 
subsamples have significant differences in sensitivity in 

37 Although recent observations suggest some (a t least small 
degre e) alignment of AGN and their host galaxies (ILagos et alj 

mm . 
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Fig. 14. — Distribution of the Compton-thickness parameter 
T = L 2 - 10k ev/L[o m ] A5007 for an X-ray+[0 m]-selected sam- 
ple of AGNs (black line). The distributions for the subsamples of 
galaxies detected at 24 fim and at 70 /an are overplotted in dark 
and light gray, respectively. The most absorbed systems (with 
lower values of T) are preferentially detected at 70/im. This sub- 
sample displays a much broader distribution of thickness param- 
eter and a larger relative number of absorbed AGN compared to 
both the 24/xm sample and the parent X-ray +[0 III] -selected sam- 
ple. Arrows at the top of the figure indicate the mean and stan- 
dard deviation for nearby AGNs that are either obscured, Type 2 
(0,59±1. 06 dex; in green) or u nobscured, Type 1 (1.59±0.48 dex; 
in black; IHeckman etWiWOS ) . 

terms of the SFR of their hosts but we otherwise con- 
strain this portion of the analysis to galaxies that have 
log(M*) > 10 10 Mo and P(AGN)>30% 38 . With compa- 
rable masses and fairly secure AGNs, the 70 /mi subsam- 
ple will include galaxies with the highest sSFRs whereas 
the 24 /mi and X+ [O in] subsamples probe down to lower 
sSFRs given a range of masses and redshifts. 

The distributions of the Compton-thickness param- 
eter (logT = log(i 2 -iofe e y/i[o in])) between the FIR 
(70 /mi), MIR (24 /itm) and X+[0 in] samples are shown 
in Figure HU The FIR-selected subsample tends to host 
more absorbed AGNs than either of the comparison sam- 
ples. If the parent population of the FIR (light gray) and 
X+[0 in] (black line without the light gray area) sub- 
samples were the same, the probability of drawing both 
distributions would be 0.0012% according to a KS test. 
It remains very low (0.038%) for the distributions of the 
FIR- and MIR-selected subsamples. The distribution of 
the latter is intermediate between the X+[0 in] subsam- 
ple and the most IR-luminous (detected at 70 /tm) sam- 
ple. Clearly, the majority of the most absorbed systems 
(with Compton-thickness log(T) < 0.3; Appendix |A| 
tend to be infrared-luminous given their high detection 
rate at 70 /im. 

The overall trend is more frequent AGN X-ray absorp- 
tion in galaxies with higher sSFRs. In order for this trend 
to be caused by geometry of the putative torus from AGN 
unified models, the small-scale obscuring torus would 
have to be physically and/or optically thicker in FIR- 
luminous galaxies with elevated sSFR. Another possibil- 
ity is that (part of) the gas reservoir available to sustain 



We set P(AGN)= 1 if the source is selected as an X-ray AGN. 



the elevated SFR also acts as an important absorber for 
the central AGN. This situation could be increasingly 
common at high er redshifts as the AGN hosts show an 
increase in SFR (IMullanev et alj |2010t IShao et aHl201Ct 
iMullanev et al.l l2012f ). Our findings imply that it is 
not just the total SFR (or Lm) that is relevant but 
rather the sSFR. This points to another factor explain- 
ing the apparent connection between AGN X-ray absorp- 
tion and host galaxies. We postulate that it could be in 
part because the sSFR can trace the gas fraction across 
the sSFR sequence ( but no longer when reaching the 
more extreme sSFRs; Mag dis et alJ 120121 ) . and because 
for the more extreme objects above the main sequence, 
the enhanced sSFR could reflect an increased gas den- 
sity. The latter can be caused by major galaxy mergers 
(e.g.. iNaravanan et al J 12001 Uuneau et alJl2009t ). which 
are common for high -redshift IR-selected ga laxies above 
the sSFR sequence (|Kartaltene et alJ l2012f ) and which 
are postul ated to result in a n increased star formation 
efficiency (Da ddi et alJl2010T ) due to the "compacting" 
of the star- forming gas into a dense central starburst. 
Recent observations showed a correlation between AGN 
X-ray ab sorption and ~ 100 pc sca l e star burst in nearby 
Seyferts (|Diamond-Stanic fe Riekdl2011[ ). which can be 
broadly consistent with either of the above scenarios 
(torus- host galaxy link, or obscuration by ~ 100 pc scale 
gas). Some simulations pred ict a multi-scale connection 
where the black hole fue ling (|Hopkins fe Quatacrt l2010T) 
and/or torus properties (jHopkins et alJl2012h are related 
to the nuclear SFR via gas inflow albeit in a complicated 
and perhaps stochastic way. These possibilities are both 
very interesting as they differ from the typical AGN uni- 
fied models where the very central region is assumed to 
be somewhat decoupled from its large-scale surround- 
ings. 

Complementary findings suggest that host-scale obscu- 
ration may hinder the detection of optical l ines in hard 
X-ray selected but optically-dull AGNs (e.g.. lRigbv et alJ 
2006). These tend to reside in disky edge-on hosts and, 
if all optical signatures are lost, would be missed with 
the MEx diagnostic diagram. The host obscuration that 
we find here must have a smaller covering fraction of the 
narrow line regions in order to detect optical nebular line 
signatures in many cases. 

6.3. Properties of Absorbed and X-unabsorbed AGNs: 
An Evolutionary Sequence? 

The stellar mass distributions of the weak and ab- 
sorbed AGN samples are similar to one another while 
the mass distribution of X-ray AGN is skewed toward 
higher stellar masses. According to Figure Eta), the 
secure X-ray AGN all have a high stellar mass M* > 
10 10 6 Mq, and seem to be more massive on average 
than the global AGN population. This suggests that X- 
ray selected AGN will tend to be more easily identified 
in more massive hosts. However, the statistics are poor 
with only 12 X-r ay AGN on the MEx diagram. Using 
a larger sample, IMullanev et alJ (|2012t ) found that X- 
ray selected AGNs always lie in high stellar mass hosts 
((-M*) = io 10 - 7 - 10 - 8 Mq) across a broad range of redshift 
(up to z ~ 3). This may be due to the presence of more 
massive black holes in these systems although the au- 
thors only found a weak correlation between stellar mass 
and X-ray luminosity over the 10 42 — 10 44 erg s -1 range 
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(M* oc L { x ,7) ). lAird et al.l (|2012h explain this feature 
with observational bias such that more massive black 
holes are observable at a broader range of Eddington 
ratios than less massive ones, which presumably reside 
in less massive hosts. While this selection bias will be 
especially true for X-ray selection, it is less severe with 
optical AGN diagnostics, which are sensitive down to 
much lower accretion rates. This may explain why the 
MEx method uncovers AGN down to lower stellar masses 
than the less sensitive X-ray method. However, every 
AGN selection is incomplete. Thus, it is of great interest 
to combine multi-wavelength indicators for the sake of 
completeness. 

Relative to X-unabsorbed AGNs, the absorbed AGNs 
have similar bolometric luminosities but somewhat 
smaller stellar masses (and possibly smaller mass black 
holes). The absorbed AGN hosts are thus growing both 
their BHs and their stellar content faster (as they tend 
to have higher sSFR), in agreement with the results from 
iTanakal (|2012f ). This is unlikely to result from an arte- 
fact due to [O in] being strongly contaminated by star 
formation because the [O in]-to-IR ratio is elevated in a 
similar fashion in both the absorbed and X-unabsorbed 
AGN hosts (right panel of Figure [5]) . 

Those systematic differences between the absorbed and 
unabsorbed AGN hosts may correspond to an evolution- 
ary sequence where absorbed AGNs tend to occur at the 
beginning of a growth cycle with higher specific growth 
rate and elusive X-ray signatures due to gaseous sur- 
roundings (on possibly a range of physical scales). X- 
ray unabsorbed AGNs would then correspond to more 
mature systems (consistent with higher stellar masses 
of the hosts) after some of the gas responsible for fu- 
eling and/or obscuring the central region has been con- 
sumed, pushed aside or removed. This idea is similar 
to th e well-known major merger scenario (| Sanders et al.l 
119881 described in more details in the next section), but 
could be more generic and also apply to AGNs with in- 
termediate luminosities in isolated but initially gas-rich 
galaxies, assuming that even in those isolated systems, 
the obscuring material could be gradually removed by 
the AGN and/or star formation. The less extreme AGN 
luminosities compared to the quasar regime mean that 
major galaxy mergers may not be required in all cases. 

6.4. Physical Origin of High AGN Fractions 

This study finds a larger AGN fraction (37%) among 
IR-luminous, star-forming galaxies than previous stud- 
ies. Com bining multiwave l ength d iagnostics from X-ray 
to radio, IKartaltepe et al.l (|2010af) found that the AGN 
fraction rises steeply with IR luminosity and / or redshift 
(see the filled black circles on the left panel of Figure E]), 
but their values are lower than ours with an overall frac- 
tion ~ 10 — 20% over the redshift and luminosity ranges 
of int erest in this paper. Similarly, ISvmeonidis et ahl 
(2010) studied the occurrence of AGN in 70/Ltm galaxies 
over a similar redshift range and found an AGN frac- 
tion ~ 13%. The main difference between these previ- 
ous investigations and the current work is the use of the 
MEx diagnostic diagram, which unveils the majority of 
the AGN including absorbed and weak systems. Some 
absorbed AGN were successfully selected by previous au- 
thors with, e.g., mid-IR color diagnostics, but the latter 
are not as sensitive to the presence of active nuclei as the 



MEx method and are therefore less complete. 

What are the physical causes of such a frequent occur- 
rence of AGN? Is there a connection between AGN trig- 
gering and the elevated SFR in their host galaxies? Or 
between AGN absorption and their host galaxies? The 
results from this work imply that active SMBH growth 
occurs in parallel with active star formation. The ele- 
vated SFRs are probed by the FIR emission and there is 
a connection between the infrared luminosity and the oc- 
currence of AGN. One interpretation that fits all of these 
observational trend s is the major merger scenario (e.g., 
I Sanders et al.lll988l ). Major mergers of gas-rich galaxies 
are expected to fuel both starbursts and AGN, with the 
AGN originally deeply buried in large amounts of gas be- 
fore emerging as X-ray detected AGN. Qualitatively, our 
observations fit a merger scenario with simultaneous star 
formation and AGN activity. In contrast, several authors 
have pointed out that the morphologies of X-ray AGNs 
do not ap pear to be mergers more often than non-X- 
ray AGNs (ICisternas et aLlboilt ISchawinski et al.ll201lt 
iMullanev et al.l l2012t iKocevski et al.l 120111 )" However , 
these studies did not include the X-ray absorbed systems 
identified by the MEx method: it is possible that the X- 
ray selected AGNs are not associated with mergers (or 
that they represent a different evolutionary stage) while 
the more absorbed AGNs are more directly triggered by 
mergers. 

An interesting feature from this work is that the rise 
of the global AGN fraction with Lm behaves similarly 
to the rise in the fraction of galaxies undergoing mergers 
and interactions and/or having an irregular morphology 
at z ~ 1 (dashed and dott ed lines on Figure HI taken 
from IKartaltepe et alJ l20 1 2) 39 . In contrast, the fraction 
of X-ray unabsorbed AGNs is much lower and does not 
track well with the fraction of mergers and/or irregular 
galaxies. While somewhat indirect, this argument sug- 
gests that X-ray selection of AGN may miss the connec- 
tion between AGN occurrence and galaxy interactions 
because of a possible mismatch between the timescale 
of AGN X-ray observability of the visibility timescale of 
mergers. Future studies of the connection between AGNs 
and their host galaxies should include both absorbed and 
unabsorbed AGNs to provide us with a more conclusive 
test. 

Despite the apparent consistency of a merger picture, 
major galaxy mergers may not be required to explain 
our results. For instance, an alternative way to obtain 
a high incidence of AGNs and also a high incidence of 
AGN absorption could be large-scale disk instabilities in 
gas-rich star-forming gal axies, as exemplified b y high- 
redshift clumpy galaxies (jElmegreen et al.ll2007l ). These 
systems are more common at highe r redshift (z > 2) 
and may be explained by simulations (jAgertz et al.ll2009l) 
in the conte xt of cold stream accretion cosmological 
models (e.g.. iDekel et al.l l2009f ) . Clumpy galaxies per- 
sist dow n to intermediate redsh ift in intermediate mass 
galaxies (jElmegreen et al.l I2009D . Unstable (or clu mpy) 
disk gala xies are predicted dBournaud et al.l l2011h and 
observed (jBournaud et al.ll2012l ) to efficiently fuel AGNs, 

39 These fractions were calculated in a slightly different redshift 
range (0.8 < z < 1.2 rather than 0. 3 < z < 1) but are simil ar to 
what is expected at lower redshifts (Kartaltcpc et al. 2010b) and 
so should not be significantly different at the exact redshift range 
of interest here. 
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including a high likelihood of absorption of X-ray signa- 
tures. 

Irrespectively of the morpholo g y and kinematics of the 
host galaxies, work by [Tanaka (2012) on weak AGNs 
suggests that, like stronger AGNs, they occur predom- 
inantly in massive hosts (> 10 10 M Q ). Those authors 
argued that the stellar mass of the host galaxies appears 
to act as a "switch" in triggering AGNs, but the under- 
lying mechanism remains unknown given that they do 
not find dependencies on galaxy morphology. The weak 
AGNs from this Paper would also be consistent with a 
"switch" at ~ 10 10 M but we do not probe to sufficiently 
low masses to really characterize that feature. Further- 
more, the MEx diagram was designed with SDSS priors 
and so will only be sensitive to the stellar mass range 
overwhich the bulk of low-redshift SDSS AGNs can be 
identified with the traditional BPT diagrams. 

While we cannot distinguish between major galaxy 
mergers, large-scale violent disk instabilities or other 
triggering mechanism without additional observations, it 
is clear that a complete picture of the growth of galaxies 
and their central SMBHs must include a gas-rich phase 
with concurrent SF and AGN, where the central AGN 
is often absorbed even for moderate-luminosity AGN 
(10 43 - 4 <Z b ol(AGN)< 10 45 ergs- 1 , i.e., below the quasar 
regime) . 

Lastly, while the AGN fraction of the intermediate- 
redshift FIR-selected sample is similar to that in nearby 
FIR-selected galaxies, the number of IR-luminous galax- 
ies is much larger at higher redshift. This implies that the 
global AGN fraction in star-forming galaxies was higher 
in the past. A greater availability of gas at earlier epochs 
may be responsible for the more common occurrence of 
highly star-forming (thus IR-bright) galaxies. The high 
fraction of AGN, even if consistent with that at low- 
redshift, implies that a much larger population of AGN 
- and especially X-ray absorbed AGN - in IR-luminous 
hosts exists at higher redshift. We may be starting to un- 
cover the missing AGN population that h as been inferred 
from c osmic X-ray back g round studies ([ Comastr i et al.l 
19951 iMushotzkv et all 120001: [Alexander et al.l 120031: 
Bauer et alJl2004t iGilli et al.ll2007ft . 

7. SUMMARY 

In this paper, we investigated the incidence of AGN 
among star-forming galaxies at intermediate redshift 
(0.3 < z < 1). The AGN fraction and AGN char- 
acteristics (X-ray unabsorbed, X-ray absorbed, intrin- 
sically weak) are function of host galaxy properties. 
AGNs w ere identified based on four diagnostics: MEx 
diagram (1 Juneau et al.l l20lTb . X-ray criteria (similar to 
iBauer et a l. 2004), IRAC color-color (jStern et al.ll2005h 
or IRAC [3.6]-[8.0] color versus redshift (Section 
The last two are similar but not identical. Our main 
findings are: 

1. Combining all AGN diagnostics, the global AGN 
fraction is 37(30) ±3% in the 70^m selected 
galaxy population including (excluding) AGNs 
less luminous than Lx = 10 42 ergs -1 (£& 2 = 
10 43 4 ergs -1 ). This AGN fraction is around a fac- 
tor of two greater than previous results in similar 
infrared luminosity and redshift ranges. 

2. The fraction of star-forming galaxies hosting an 



AGN increases as a function of Lir. The AGN 
fractions are very similar to those in nearby (0.05 < 
z < 0.1) FIR-selected galaxies suggesting mild or 
no evolution since z = 1. 

3. The differences between the higher fraction pre- 
sented here relative to previous studies at compa- 
rable redshifts may be accounted for by (i) heav- 
ily absorbed AGNs, and (ii) intrinsically weaker 
AGNs. Absorbed and/or weak systems are more 
difficult to detect but can be identified thanks to 
the high sensitivity of the MEx diagnostic. 

4. The fraction of galaxies hosting an AGN appears 
to be independent of the specific star formation 
rate of the host, and remains elevated both on the 
sSFR sequence and above. In contrast, the frac- 
tion of AGNs that are X-ray absorbed increases 
substantially with increasing sSFR, possibly due 
to an increased gas fraction or gas density of 
the main/central star- forming regions in the host 
galaxies. 

5. AGNs with the most X-ray absorption (Njj > 
10 24 cm~ 2 ; inferred from log(£x/£[o m]) < 0-3) 
tend to reside in IR-luminous galaxies (detected at 
70 fim). Adding less IR-luminous AGN host sam- 
ples (undetected at 70 /im) mostly contributes un- 
absorbed and mildly absorbed AGNs (Section l6.2|) . 
Together with item 4, this result suggests a con- 
nection between the host galaxy's gas content and 
absorption of the central engine, which can be 
achieved either through absorption of the AGN X- 
rays by the ISM of the host galaxy or by a physical, 
multi-scale link between the host galaxy's gas con- 
tent and the obscuring torus. 

6. AGN radiation is likely responsible for the high- 
est 24-to-70 /xm flux ratios in the most extreme 
cases, which tend to be X-ray and optically setected 
AGNs. However, other factors must be at work to 
explain the spread of the mid-to-far IR color among 
star-forming galaxies and galaxies with a weak or 
absorbed AGN. 
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APPENDIX 
AGN X-RAY ABSORPTION 

This Paper uses the Compton-thickness parameter as a proxy for AGN X-ray absorption (Sections l5.2l and l5.5[) . While 
this approach may not yield accurate measurements of the absorbing column density (Nh), it allows us to separate 
the most heavily absorbed cases from the rest of the AGN pop ulation (X-ray unabs orbed and mildly absorbed) . We 
demonstrate this with the analysis of nearby AGNs presented bv lBassani et a l. (1999, hereafter B99). The B99 sample 
is not statistically complete but it spans a broad range of X-ray absorption properties and includes measurements of 
[O ill] A5007 and 2 — lOkcV fluxes, iron Ka equivalent widths, and Nh from fitting the X-ray observations with a 
photoelectric cutoff model (assuming no components from reflection or stellar emission). In some cases, the authors 
note Compton-thick signatures with Nh > 10 25 cm~ 2 if there is also absorption of X-rays at energies above lOkeV, or 
with Nh > 10 24 cm~ 2 if there is no data available beyond lOkeV to rule out an absorbing column between 10 24 cm~ 2 
and 10 25 cm~ 2 (see th eir T able 2). 

As shown in Figure |AT| the hard X-ray (2-10 keV) to [O III] A5007 luminosity ratio drops when Compton-thickness 
is reached at Nh ~ 10 24 cm~ 2 (filled black circles). However, the most heavily-absorbed systems can be mistaken 
for low-column density cases with simple X-ray spectra models that do not account for additional signatures such as 
the equivalent width of the Fe Ka emission line (which becomes high in cases of heavy absorption) or the absorption 
of even higher energy photons (>10keV). This feature is obvious in Figure [ATI where the two estimates of Nh for 
Compton-thick AGNs (see B99 for details on Nh calculations) can differ by 5 dex and typically by 3 dex. This shows 
the failure of the photoelectric cutoff estimate of Nh in cases of the heaviest absorption. 

The Compton-thickness parameter tends to be low for the most absorbed systems. We adopted the threshold 
log(Lx/Ltp ln ])<0.3dex to classify two AGNs that met the X-ray AGN criteria as X-absorbed AGNs rather than 
X-unabsorbed AGNs fSection l5.2l and Table 1). However, it is possible that other AGNs in the X-unabsorbed category 
have in fact significant absorption as we do not use X-ray spectral signatures due to the low number of X-ray counts 
(and associated uncertainties) for AGNs in our sample. In that sense, the absorbed fraction is a conservative lower 
limit and could be even higher if we moved more AGNs from the X-unabsorbed to the X-absorbed category. 

In order to facilitate the comparison with the sample of interest in this work, we used the [O in] A5007 fluxes 
uncorrected for dust obscuration. Applying the correction from the Balmer decrement method, as done by B99, 
would increase the [O in] A5007 fluxes by 1.6 (1.0) dex on average (median) with some cases presenting more extreme 
obscuration (up to 2.7 dex). In particular, two of the three Compton-thick AGNs with \og(Lx / L[ ln ])>0.3dex have 
extreme Balmer decrement values (>10), which implies that the observed [O in] A5007 line is underluminous and leads 
to a higher Compton-thickness parameter. 
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